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ABSTRACT 

The virtual mass of circular cylinders immersed in a uniformly 
stratified nonhomogeneous fluid medium has been investigated both theo- 
retically and experimentally. The experimental equipment consisted of 
a reservoir filled with one inch thick layers of salt water solution. 

The density uniformly decreased with elevation (stable stratification). 

The fluid motion about the cylinder was provided by the standing waves 
created with a plunger. The resulting velocities and accelerations 
about the cylinder were harmonic. The theoretical analysis consisted 
of the calculation of the inertia coefficient through the evaluation of 
the total kinetic energy of £ fluid jmedium. Use was made of the psuedo- 
potential flow solution for the stratified flow about a circular cylinder. 
The results have shown that the inertia coefficient for circular cylinders 
immersed in a uniformly stratified media is about 25% larger than that for 
a cylinder immersed in an homogeneous medium. 
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DEFINITION 

2 2 
Area, ft or in 

Amplitude of wave, inches 

Inertia coefficient 

Stream function coefficient 

Stream function coefficient 

Force, pounds 

2 

Acceleration of gravity, ft. /sec 

Mean height of free surface of water above bottom 
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Bernoulli "constant 11 

Wave number, sec ^ 
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u Horizontal component of velocity in homogeneous 

medium, ft/sec 

u^ Horizontal component of velocity to left of cylinder, ft/sec 

U£ Horizontal component of velocity to right of cylinder, ft/sec 

u' Horizontal component of velocity in stratified medium, ft/sec 

V Total velocity, ft/sec 

v^ Radial velocity in polar coordinates, ft/sec 

Vq Tangential velocity in polar coordinates, ft/sec 
w Vertical component of velocity in homogeneous medium, ft/sec 

w 1 Vertical component of velocity in stratified medium, ft/sec 
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Acceleration, ft/sec 
Vertical component of vorticity 
An angle in degrees 
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Stream function 
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I. Introduction. 

The purpose of the present investigation is the determination of 
the virtual mass coefficient for circular cylinders immersed in a uni- 
formly stratified nonhomogeneous medium. Virtual mass, which is a 
combination of the mass of the body and the added mass resulting from 
the acceleration of the fluid particles, is of considerable importance 
in the determination of the resistance of bodies accelerating from rest 
or decelerating to a full stop. It is for this reason that numerous 
studies were conducted to elucidate the concept of virtual mass and to 
determine the appropriate virtual mass coefficients for various shapes 
of bodies in homogeneous fluid media. 

It is necessary to briefly describe various definitions of the 
added mass concept prior to a description of the work undertaken here- 
with. As a body is accelerated from rest in a fluid otherwise at rest, 
the fluid medium surrounding the body is also accelerated. The signif- 
icant difference between the acceleration of the body and that of the 
fluid particles is that while the solid body has constant acceleration, 
the fluid particles have accelerations varying from that of the solid 
body to zero at an infinite distance from the body. In order to deter- 
mine the acceleration of the solid body for a given force F exerted on 
the body, it is, therefore, necessary to know not only the mass of the 
body but also the integral contribution of the fluid medium surrounding 
the body. If a volume integral is taken of the local accelerations 
times the infinitesimal fluid masses and expressed in terms of a lump 
sum of fluid mass times the constant acceleration of the body, then that 
lump sum of the fluid mass may be regarded as the mass attached to the 
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solid body or as increasing the mass of the solid body by this additional 
mass of fluid. It is for this reason that the mass contributed to the 
solid body is called the added mass and the sum of the two masses is 
called the virtual or apparent mass. In other words, the virtual mass 
is the quotient of the force required to produce the accelerations 
throughout the fluid divided by the acceleration of the body. In 
mathematical terms, Newton’s second law may thus be written as 

F » (M + m)<* 

in which M is the mass of the solid body, m the added mass, and cH, is 
the acceleration of the solid body. Obviously, in the case when the 
body is immersed in a fluid medium of very small density such as air, m 
is negligibly small and thus the equation of motion may be written in 
its customary form of F = M<* . In the case where the body is immersed 
in a fluid of relatively large density such as water, then m becomes 
quite significant and its neglect may lead to extremely large errors 
depending on the actual mass of the body. As an example, it is easy 
to verify that the acceleration of an air bubble will come out to be 
infinitly large if the added mass is neglected, whereas its true value 
is twice the acceleration of gravity as would be found through the in- 
corporation of the added mass into the equation of motion. 

In operations in large liquid bodies such as in the oceans, one 
must take into consideration not only the contribution of the added mass 
in determining the forces acting on the body, but also the effect of the 
variations in density of the fluid medium on the added mass itself. 

The studies carried out on added mass in the past were confined to homo- 
geneous fluid media. The advent of underwater technology and the recent 
progress made in understanding the motion of stratified fluids made it 
clear that the engineer must have a better understanding of the virtual 
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mass of various shapes of bodies in nonhomogeneous fluids. It is 
apparent that there are an infinitely large number of combinations of 
shapes of bodies and types of liquid nonhomogeneities to be considered. 
For a scientific approach to the problem, however, one must consider 
the experimentally and mathematically amenable cases such as the motion 
of a right circular cylinder in a uniformly stratified fluid medium. 

It is for this reason that the project described herein investigates the 
virtual mass of three sizes of circular cylinders in two different but 
uniformly stratified salt-water solutions through the use of standing 
waves. The results show, as anticipated theoretically, that the vir- 
tual mass coefficients in stratified fluids are about 25% larger than 
those in homogeneous fluids. 
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II. Previous Theoretical and Experimental Studies. 

The concept of virtual mass was probably first recognized by DuBuat 
in 1786 while studying the sources of errors in the motion of pendulums. 
Bessel, in 1826, also recognized the effect in his studies of the same 
problem. Analytical solutions to the problem were found by Poisson and 
Green in 1832 and 1836 . Froude, in his studies of ship oscillations in 
1861, also noticed the added mass effect and took it into account in his 
inertia coefficient. 

Since 1920, there have been a number of papers and texts published 
concerning added mass, virtual mass, or inertia coefficients. Lamb (1)* 
developed some of the concepts and mathematical background of the virtual 
mass effect. Landweber and Yih (7) compared the solutions of Lagally, 
Taylor, and Kirchoff showing their correlation and simplying their proofs. 
Texts by Rouse (3) and Robertson (2) give clear definitions and theoreti- 
cal developments of the phenomenon. 

Much experimental work on virtual mass has been done in the past 
15 years. Keulegan and Carpenter (5) determined inertia coefficients 
for circular cylinders and flat plates through the use of standing 
waves. McNown (9) carried out similar studies for lenticular cylinders 
and flat plates. Sarpkaya (13) conducted studies on the added mass of 
lenses and parallel plates. Brater et. al. (8) dealt with the determina- 
tion of the virtual mass coefficients of barges of different geometries. 
All of these studies have tended to support the theoretically determined 
virtual mass coefficients for low velocity and small amplitude motions. 

^Numbers in parenthesis refer to bibliography Section VIII. 
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As the velocity and amplitudes were increased, departures from the 
theoretical values were observed. The theoretical analyses invaribly 
assumed irrotational flow in an ideal, homogeneous, incompressible 
fluid. Thus it was to be expected that experimental values would 
differ from those obtained theoretically as separation, vortex forma- 
tion, and turbulence accompanied the flow following the initial periods 
of motion. Sarpkaya (13,14) carried out extensive theoretical and ex- 
perimental studies of the virtual mass and time-dependent resistance 
for flows with vortex separation. 

The previous experimental studies carried out with homogeneous in- 
compressible fluids with bodies either undergoing oscillations of small 
amplitude or accelerating from rest have shown that the theories of 
Lagally, Taylor, and Kirchoff are substantially correct provided that 
the fluid motion does not continue long enough in one direction to 
cause separation. 

The work on the corresponding problem in non-homogeneous fluids, 
however, has not been as prevalent. Yih (4,12) and Long (10,11) ap- 
pear to have done most of the work in this field in the past ten years. 
The most important discovery with regard to virtual mass effect is 
the phenomenon of ,r blocking" discussed by Yih (4). 

The experiments of Long (11) and Yih (12) indicated that for 
small amplitude motion stagnation zones occur upstream and sometimes 
even downstream from a barrier or a moving body. This phenomenon of 
stratified flow would seem to manifest itself in the form of a "frozen" 
volume of fluid following the motion of a solid body. Theoretically, 
there appears to occur a critical regime in which blocking takes place. 
If the stratification is slight, or the motion rapid, the effect is 
not noticed. In the present study, the blocking effect per.se. was 
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not studied, however the relationship between this effect and the virtual 
mass effect was necessarily taken into account. 

From the foregoing, it should be apparent that the blocking would 
increase the virtual mass coefficient. The validity of this argument 
has been shown through the experiments presented herein. 
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III. Experimental Equipment and Procedures, 

Equipment 

For the experimental study, a one inch thick plexiglaa tank with 
outside dimensions 24" by 24" by 36" was constructed. A inch hole 
was drilled and tapped for pipe threads in the center rear of the tank 
bottom. A three inch diameter circular aluminum flat plate was install' 
ed with a 3/16" clearance over this hole. The purpose of the hole was 
to permit slow filling of the salt-water in layers of varying density 
from the tank bottom. 

A pressure tap was drilled in the side of the tank 4 inches from 
the top and centered 12 inches from either side of the tank. This tap 
was made for use with a differential pressure transducer in recording of 
wave heights. A mounting bracket for the transducer was installed on 
the outside of the tank in a position such that the aero level of the 
transducer was located at the level of the pressure tap. One side of 
the transducer was connected to the tank tap and to a length of tygon 
tubing in the shape of a standpipe. The standpipe was used as a cali- 
bration device. The other side of the differential transducer was left 
open to atmosphere for the purpose of recording gage pressures, 

A heavy steel box frame was constructed around the tank and its 
supporting table. The base of the frame rested on the concrete floor of 
the building. The top of the frame was securely bolted to a concrete 
pillar. The frame was designed so thst the height of the cylinder arm 
could be adjusted in two inch Increments, This was accomplished by means 
of an aluminum bracket which bolted to both sides of the steel frame, 

The frame was drilled and tapped at two inch intervals from the level of 
the top of the tank to the top of the frame, (See Fig, 2), 
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Figure 2. A schematic drawing of the test apparatus 
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The cylinder arm was attached to the bracket by a gimbal bearing 
with only one degree of freedom. The upper extension of the arm was 
clamped to the steel frame with a "C" clamp. The cylinders themselves 
were attached to the arm by another aluminum frame. This frame consist- 
ed of two 1" by 1/8" streamlined aluminum arms, 20" long. The arms were 
bolted to a top piece which consisted of 1" by V* aluminum bars bolted 
together in the form of a channel beam. The top piece was 21. 5" long 
as were the cylinders. Thus, when the cylinders were bolted to the 
arms, a 20" by 21. 5" rectangular rigid form was produced. This was in 
turn bolted to the cylinder arm at the center of the top piece. 

The cylinder arm was a 3/4" by 7/8" steel rod. A one inch test sec- 
tion was notched into the arm at a point six inches above the attach- 
ment point of the cylinder frame. The cross section of the test section 
was 3/8" by 7/8". A four gage bridge was formed at this section using 
Microsensor strain gages with a gage factor of 150. The bridge was 
completed at the test section and the leads brought out in a shielded 
cable to one channel of a two channel Sanborn amplifier/recorder system. 

Calibration of the bridge was accomplished through a pulley, weight, 
and string system. One pulley was attached to a heavy brass block and 
set into the empty tank. Another pulley was clamped to the aluminum 
cross-bracket. A string was attached to the middle of the cylinder 
and passed over both pulleys. A weight pan was built to attach cali- 
bration weights to the string. Thus a horizontal force was produced on 
the cylinder equal in magnitude to the weight of the pan and weights. 

(See Fig. 3). 

The tank itself was mounted on a heavy 4’ by 4* table with a hole in 
the center rear of the table. A six inch long and %-inch diameter pipe 
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was attached to the hole in the bottom of the tank, and supported by a 
pipe bracket bolted to the table. This pipe extended below the table 
and was fitted with a valve and tubing fittings for filling or draining 
the tank. (See Fig. 2). This valve is shown as valve A. The filling 
tubing ran to valve B located at the front of the table. Tubing from 
B led to a tee connection at D which was the supply spigot from the 
building water main. From spigot D tubing was led through another tee 
either to the filling bucket or to a washdown line. Valves C and E were 
isolation valves for their respective lines. 

A simple filling and mixing bucket was made by punching a hole in 
the bottom of an old five gallon kerosene can which had been thoroughly 
cleaned. A fitting was attached to the bottom of the bucket and tubing 
run from there to the filling system as shown on Fig. 2. Thus, as can 
be seen from the diagram, by properly manipulating the valve arrange- 
ments, the tank could be filled either from the building supply with 
fresh water, or from the bucket with salt-water. 

Salt-water solutions were made in the bucket by filling the bucket 
with fresh water to a height corresponding to a one inch layer in the 
tank, then adding measured quantities of salt. In order to increase the 
pressure head while filling from the bucket, the bucket was raised about 
six feet above the level of the table. The salt solution could then 
drain by gravity into the bottom of the tank. 

A standing wave was created in the tank by means of a 20 inch by 8 
inch piece of % inch plywood. Using this as a paddle of sorts, a wave 
was gently started by hand. The amplitude of the wave was increased by 
applying a slight lift or push to the paddle in phase with the natural 
motion of the standing wave in the tank. The plywood was left floating 
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Figure 3 Calibration System 




in the tank to help damp out any secondary waves, and because lifting the 
paddle out of the tank tended to create secondary waves disturbing the 
basic standing wave pattern. 

All readings were taken on a two channel Sanborn amplifier/recorder. 
This instrument consisted of a Hewlett/Packard model 7712-B recorder 
with Sanborn model 350-1100C preamplifiers. 

Procedures 

Prior to conducting the experiments a careful calibration procedure 
was followed. The first item requiring calibration while the tank was 
empty was the strain gage bridge and recorder system. This calibration 
was accomplished by the pulley system previously described and shown in 
figure 3. With no weights or weight pan on the string, and all tension in 
the string relaxed, a careful zero balance was obtained on the recorder. 
Then an accurate weight of 0.1 pounds, including the weight of the pan, 
was suspended from the string. The needle deflection on the recorder 
was adjusted by means of the preamplifier gain control so that at an 
attenuation setting of 10, a deflection of 25 millimeters (half of full 
scale) was indicated. The weight was then removed and the system zero- 
balanced again. The above procedure was followed again and checked to 
insure that the deflection for the same weight and attenuation setting 
remained constant. All higher attenuation settings were then checked to 
insure linearity. When this had been accomplished for 0.1 pounds, the 
weight was changed to 0.2 pounds and the procedure followed again. Once 
the gain was set during the initial calibration, it was locked in place. 
Calibration checks were made on both sides of the cylinder, that is the 
pulley system was shifted to the other end of the tank and calibration 
repeated for force in that direction. Calibration remained the same for 
either side of the cylinder. 

25 



The tank was then filled to a depth of 19 inches with fresh water. 
This allowed for one inch of water above the pressure tap in the side 
of the tank. The transducer pressure side was filled with water from 
the tank and all air pockets removed from the tubing and transducer 
cavity. The transducer output was connected to the other channel of 
the recorder. The recorder was again zero balanced with both the tank 
valve and the standpipe valve open. The tank stop was then shut and a 
water column of one inch above tank level was impressed on the trans- 
ducer by means of the tygon tubing standpipe. The amplifier gain was 
adjusted for half of full scale deflection with the attenuation setting 
at 2. The gain on that channel was locked and checks of linearity were 
made similar to those conducted for the strain-gage arm. Water heights 
of up to one inch above and below the level of the tank were tested, 
and the calibration was found to be linear. 

Following the completion of the calibrations, a complete system 
check of the experimental apparatus was made. With fresh water in the 
tank, runs were made with the 1.0, 1.5, and 2.0 inch diameter cylinders 
to see how closely the experimentally determined homogeneous fluid 
virtual mass coefficients agreed with those predicted theoretically. 

To do this, a standing wave was generated as previously described. Once 
the fundamental standing wave with no harmonics had been set up, the force 
acting on the cylinder and the wave heights at the end of the tank were 
simultaneously recorded. 

Immediately after each experiment, the transducer calibration was 
checked. It was found that the gain setting needed no changes through- 
out the series of experiments. The water was then drained out of the 
tank and the strain-gage am calibration checked again as above. Again, 
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no change was required from the initial calibration setting in spite of 
changing of the cylinders. 

When all the runs had been completed with fresh water, the water was 
drained out of the reservoir, calibrations were repeated, and the strati- 
fied layers were introduced into the tank. This was accomplished by 
slowly filling the water-salt solution from the bottom of the tank through 
the filling system previously described. The first layer to be intro- 
duced was two inch layer of fresh water. The solutions for the second 
and succeeding layers were prepared by filling the bucket and adding a 
measured quantity of salt (10 milliliters or 30 milliliters) and mixing 
thoroughly. The quantity of salt was Increased by equal amounts for each 
succeeding layer, that is, 10 ml. for the second layer, 20 ml. for the 
third, 30 ml. for the fourth, etc. As each solution was mixed in the buck- 
et it was dyed with red, blue, or black ink, or on occasions was left 
clear for contrast. The purpose of the dye was merely to make the layers 
visible. Once the salt, water, and dye had been mixed and the salt com- 
pletely dissolved, the specific gravity of the solution was determined 
with a hydrometer. The bucket was then raised to an elevation 6 feet a- 
bove the level of the reservoir, and the valves opened to allow draining 
of the bucket into the bottom of the tank. 

The flat plate over the filling hole allowed the slow radial diffu- 
sion of the entering solution, preventing premature or undesirable mix- 
ing of the successive layers. As the new, more dense layer entered the 
bottom of the tank and spread radially from the filling hole, it lifted 
the less dense layers above it. This process was continued until the 
tank had been filled to a water depth of approximately 19 inches. Once 
the filling process was completed, the tank was allowed to sit for a day 
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Figure 4(a) Stratified Layers 




Figure 4(b) Wave Form 



r ) Q 



or so prior to the start of the tests. This waiting period allowed mole- 
cular diffusion to take place producing a uniform density gradient rather 
than a step layer effect. Each one inch layer took approximately 20 
minutes to prepare, and an hour to drain into the tank. The entire fill- 
ing procedure took a period of three days. 

Once the tank was filled with the stratified salt-water, the cali- 
bration of the transducer was checked as previously described. A stand- 
ing wave was set up in the tank. The forces acting on the cylinder and 
the wave amplitudes were recorded. Two density gradients and three 
different cylinder diameters were used in the experiments. From these 
data the virtual mass coefficients of cylinders in both the homogeneous 
and stratified fluid media were calculated. 
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IV. Theoretical Considerations. 



1. Added Mass and Inertia Coefficients (Homogeneous Fluid) 

The force exerted on a circular cylinder by an inviscid fluid 
of velocity U and acceleration dU/dt may be determined as follows: 

Letting U be the undisturbed velocity of the ambient flow, 
the velocity components for the flow about the cylinder are given by 
(See Lamb) , ,, z 



U = U (l + “S. ca*>29) 



W = - U“4 s,n Z9 



( 1 ) 



The pressure at any point in the fluid medium reduces to: 

f = 3T ( r+J £) sln ® 



( 2 ) 



Writing the equation of momentum in its familiar form 

r “ *c 

p a - J d A -v- ^ (u*-u£)dX + j ( pi “ Pa) d * (3) 

A — co -os 

and noting that for the case under consideration, » u^ at infinity, 
i.e., the flow is symmetrical, we have 




It is apparent that the force acting on the cylinder in the 
absence of a pressure gradient is given by the first term on the right 
hand side of equation (4). This corresponds to the case where the fluid 
is at rest and the cylinder is moving with the time dependent velocity U, 
Because of this obvious difference between the acceleration of flow about 
a cylinder at rest and the acceleration of the cylinder in a fluid at 
rest, two types of force coefficients are employed in the literature. 
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The normalized force coefficient for the case corresponding to the ac- 
celeration of the fluid is called the inertia coefficient, and the one 
corresponding to the acceleration of the cylinder, the added mass co- 
efficient. As will be seen later, the two coefficients differ only and 
always by a constant equal to 1.0. 

Inserting equations (1) and (2) into equation (4) and performing the 
indicated integrations, we have 

P = TT P r 0 2 iy 4- Tr P To 2 - dU (5) 

J dt J dt 

Normalizing equation (5) by "displaced mass times acceleration", 
yields 



F 



ir f 



0 dt 




( 6 ) 



where Cm is known as the inertia coefficient. It is obvious that 

in the case of a circular cylinder accelerating in a fluid otherwise 

at rest k* « 1, k* being the added mass coefficient. For other shapes 

of bodies, k* differs from unity and C reduces to: C = k* + 1. In 

m m 

passing it should be noted that for a sphere C = 1.50, provided that 

m 

the flow about the sphere is unseparated. 

In the present investigation, the cylinder is at rest and the theo 

retical value of the inertia coefficient for an homogeneous working 

fluid is C “ 2.0. This will later be compared with those obtained 
m 

experimentally and it will be shown that there is excellent agreement 
between the theoretical prediction and the measured values. 



2. Added Mass Coefficients (Nonhomogeneous Fluid). 

The equations of motion for a nonhomogeneous fluid will first 
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be summarized prior to the evaluation of the added mass coefficient for 
a circular cylinder immersed in a uniformly stratified fluid medium. 



The equation of continuity may be written as 



-f ^ W 

b X b 2 



O 



where 




( 7 ) 



( 8 ) 



The stream function 



u 






is given by 



and 




( 9 ) 



The equations of motion, in terms of u* and w 1 , become 



and 




( 10 ) 




jp- W 





(U) 



Writing the Bernoulli constant H ^ as 

H b = p + + t3 t <»> 

and the £, “ component of vortlcity f\* an 



n / w ' - 

1 bZ b X 




equations (10) and (11) may be written as 



( 13 ) 
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The foregoing equation was derived by Long (11) and Yih (4). 

Equation (14) may be used in the study of a class of psuedo- 
potential flows. One such flow case is the motion of a circular cylinder 
in a uniformly stratified fluid. The stream function for the case under 
consideration has been given by Yih (4) as follows: 



in which r^ * r/r^, r o t * ie radius of the cylinder, * U, the 

velocity of the cylinder, and 



of equation (15) by evaluating the velocity components and performing 
the integration for the kinetic energy as 




(15) 




(16) 



The added mass coefficient may be calculated through the use 




9 



(17) 



and noting that in 




(18) 



m is the added mass sought. 

The velocity components in polar coordinates are obtained from 



equation (15) as 



V r •= -r 



r ae 



Oo r * 

8 



(Si 

l r z 



cos 9 - J2l cos 30 






Cos 9 



V - -Vt' - C° r ° ~ 
“ 8 



(±s; n a-<g<,\ n 3j f c 4 Ji sm e 



Evaluating the total velocity, we have 

V a -* 1 <V = ' 

- sin 0 sin 3©} - 

— JjSL CosSccS 3&) ■+• 



5 Ik cos 9 ‘OS 30 

A. ^ ia 

C.Q,r. r /lf c „*«9 
Tpr 1 r» 

Cot-^ V 0 f fji 5,0^0 

IT 1- v rl 



JL Sin © sin 30 



5 ** 

+ C 4 f* 



(19) 

(20) 



( 21 ) 



Obviously, the integration of the terms Involving sin 9 sin 

36, cos 0 cos 30 through the use of equation (17) between the limits 

indicated shows that these terms have no contribution. Furthermore, the 

2 2 

integration of the terms involving cos 0 and sin 0 cancel each other 
out due to the opposite signs of their coefficients. The remaining terms, 
i.e., the first, second, and the last term on the right hand side of 
equation (21) yield 



34 



( 22 ) 






Hence the added mass YT\ reduces to 



ir a r ^ 4- f° 0 ^ Cs. 
^ V ?° r ° 43“ | i* 



U 



(23) 



Inserting the value of L 0 from equation (16) into equation (23) and 

Z 

dividing by the displaced mass of the cylinder i.e., IT ,we 

have the added mass coefficient k*, as 



k*= &) 



(24) 



It should be noted that for a density stratification where 
^ is decreasing with elevation (stable stratification), ^ ^/d Z 
is negative. 

The inertia coefficient Cm > which was given by Cm" k* + 1> 
now reduces to 



Lm ^ r Z \ 1 Po d z / 



(25) 



The results obtained from equation (25) will be compared in the next 
section with those obtained experimentally. It suffices to note that 
in stratified flows the inertia coefficients depend, unlike homo- 
geneous fluids, on the density of fluid, the radius of the cylinder, 
and of course on the density gradient. 
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3. 



Characteristics of the Experimental Waves. 



The region under the nodal area of a standing wave that may be 
realized in a rectangular tank provides a velocity field of simple 
harmonic motion in the velocity component U. 

Taking the x-axis in the plane surface of the undisturbed water, 
the z axis vertical and upwards, and the origin at one end of the tank, 
the surface elevation, measured from the undisturbed level, is given 
by Miche (6) to a second order of approximation in the amplitude as: 



the period of oscillation; H the depth of water; and a the semi -wave 
height, that is the mean value of the extreme end deflections in a cycle. 
The period is related to k and H by the expression 



The expressions for the particle velocities, within the order of 
the approximations considered, are given by Miche (6) as 





cosh 2kH 



where 



and 




where k = 1_ , L being the length of the basin; O' s ^' n V / T > T being 




(27) 



The surface elevation at x = 0, at any time, becomes 




(28) 
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u - 



qak cosh k(z-*-H) . i i 

■ kx cos o'! 



Cosh kH 



Cash 2k (» + H ) . , £in - t 

4 O' S'»h 4 kH sinh2kH SlriZkx S,n2fft 



" W = ^ 

O’ 



s>r>h k (2 + H ) 

cos h KH 



cos kx 



+ 1 gflh 1 Sinh 2 k (i + H) 

4 o' a;nh»kH sink 2 kH 



Cos <rt 

Cos kx sin 2c'‘t 



( 29 ) 



(30) 



At the vertical plane through the midsection of the tank, i.e., at the 
plane x » L/2 or kx ■ IT /2, the velocities become 



u 



111 coah k(»-H ) Cos a>t 
cr cosh kH 



(31) 



and 



W - - 



3 

4 






Sinh 2k (zt H) 
sinh a k H sin h kH 



sin 2 &~ t 



(32) 



Obviously, the horizontal component of velocity at the channel mid- 
section is simple harmonic. The vertical component is also simple har- 
monic except that the frequency is twice as large. The effect of the 
vertical velocity is negligibly small and decreases with both the wave 
height and the lowering of the cylinder in the reservoir. The majority 
of the experiments were made with^a^ less than 0.5 inch. For these cases 
the ratio of the amplitudes of vertical and horizontal components of 
velocity, VJm * at the cylinder location was found to be 0.004. 

The acceleration at the nodal plane may be obtained from equation 

( 31 ) 
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( 33 ) 



- Ota k cosh 

— • - j — — — tom <rt 

d t Oos k k H 



Clearly, acceleration reaches its maximum value when the velocity 
reaches its minimum value of zero and vice-versa. This observation con- 
stitutes one of the most significant points of the experimental proce- 
dure. It is a well known fact that the prediction of the forces attri- 
butable to drag and to virtual mass is quite complex for flows in which 
separation takes place. At moments of maximum acceleration and zero 
velocity, there is no wake, and the pattern of flow closely resembles 
that for irrotational flow without separation. The foregoing dis- 
cussion and analysis explain the reason for the selection of a simple- 
harmonic wave motion in the determination of virtual mass of a circular 
cylinder or, for that matter, of any other object. 



V. 



Results . 



The following tables present representative values obtained during 
the course of the experiments. The following procedure was used in the 
determination of the inertia coefficient: 




1. Inserting the measured amplitude of the wave and the appropri 
ate values of k, H, and z into equation (33), noting that the 
amplitude was taken at its maximum when cos (tf* t)* 5 1, the 
amplitude of the acceleration was calculated; 

2. Using the corresponding measured values of the amplitude of 
force, the inertia coefficient was calculated from equation 
( 6 ). 
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TABLE I 





Results of system 


check with fresh 


water 


Cylinder diameter; 1 inch 

Water depth; 19 inches 

Height of cylinder above bottom; 


4.375 inches 




RUN NUMBER ONE 








Point 


Force on 
cylinder 
(pounds) 


Amplitude 
of wave 
(inches) 


Inertia 

Coefficient 


1 


.0124 


.281 


2.154 


2 


.00976 


.25 


1.906 


3 


.0092 


.2187 


2.054 


4 


.008 


.212 


1.842 


5 


.0072 


.1875 


1.875 


6 


.006 


.141 


2.077 


7 


.006 


.141 


2.077 


8 


.0044 


.105 


2.046 


9 


.0036 


.0937 


1.876 


10 


.00264 


.0625 


2.062 


RUN NUMBER TOO 








1 


.018 


.4225 


2.080 


2 


.016 


.406 


1.924 


3 


.0145 


.360 


1.966 


4 


.0128 


.328 


1.905 


5 


.0116 


.282 


2.008 


6 


.0104 


.250 


2.031 


7 


.0092 


.235 


1.911 


8 


.0083 


.208 


1.948 


9 


.00575 


.141 


1.991 


10 


.0048 


.1158 


2.023 



40 



TABLE I (continued) 





Results of system 


check with fresh 


water 


RUN NUMBER 


THREE 








Point 




Force on 
cylinder 
(pounds) 


Amplitude 
of wave 
(inches) 


Inertia 

Coefficients 


1 




.032 


.750 


2.085 


2 




.029 


.720 


1.968 


3 




.027 


.666 


1.981 


4 




.0255 


.610 


2.043 


5 




.0235 


.579 


1.983 


6 




.0212 


.506 


2.047 


7 




.0206 


.500 


2.015 


8 




.0185 


.445 


2.032 


9 




.017 


.406 


2.046 


10 




.015 


.375 


1.955 


11 




.0136 


.344 


1.932 


12 




.0124 


.312 


1.942 


13 




.0114 


.281 


1.982 


14 




.0098 


.250 


1.916 


15 




.0094 


.219 


2.098 


16 




.0084 


.203 


2.022 


17 




.0078 


.1875 


2.033 


18 




.0065 


.1562 


2.034 


19 




.0058 


.1407 


2.015 


20 




.0052 


.125 


2.033 
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TABLE II 



Results of system check with fresh water 



Cylinder diameter 1.5 inches 
Water depth: 19 inches 

Height of cylinder above bottom: 6.5 Inches 



Point 


Force on 
cylinder 
(pounds) 


Amplitude 
of wave 
(inches) 


Inertia 

Coefficient 


1 


.0146 


.150 


1.957 


2 


.024 


.25 


1.930 


3 


.026 


.275 


1.901 


4 


.0375 


.388 


1.943 


5 


.040 


.413 


1.947 


6 


.0415 


.425 


1.963 


7 


.020 


.200 


2.010 


8 


.0218 


.225 


1.948 


9 


.0300 


.3125 


1.930 


10 


.034 


.338 


2.022 




TABLE 


III 




Cylinder diameter: 
Water depth: 18.875 
Height of cylinder 


2 Inches 
Inches 

above bottom: 8. 


5 Inches 




1 


.058 


.272 


2.097 


2 - 


.055 


.272 


1.988 


3 


.052 


.247 


2.070 


4 


.030 


.147 


2.007 


5 


.029 


.136 


2.097 


6 


.080 


.358 


2.198 


7 


.070 


.321 


2.144 


8 


.062 


.297 


2.053 


9 


.038 


.185 


2.020 


10 


.024 


.124 


1.903 
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TABLE IV 



Results in stratified salt-water 

4 

Density gradient: .000020 #/in 
Water depth: 19.5 inches 
Cylinder diameter: 1 inch 

Height of cylinder above bottom: 6.5 inches 

Point Force on Amplitude Inertia 

cylinder of wave Coefficient 

(pounds) (inches) 



1 


.0309 


.51 


2.83 


2 


.0298 


.48 


2.90 


3 


.0270 


.46 


2.74 


4 


.0258 


.425 


2.84 


5 


.024 


.40 


2.80 


6 


.0222 


.37 


2.80 


7 


.0205 


.35 


2.74 


8 


.0195 


.32 


2.85 


9 


.0172 


.30 


2.68 


10 


.0120 


.21 


2.67 


11 


.0116 


.20 


2.71 


12 


.010 


.175 


2.67 


13 


.009 


.16 


2.63 


14 


.008 


.13 


2.88 


15 


.0073 


.12 


2.84 


16 


.032 


.52 


2.87 


17 


.0289 


.495 


2.72 


18 


.026 


.43 


2.82 


19 


.0245 


.395 


2.89 


20 


.0173 


.29 


2.79 


21 


.160 


.26 


2.87 


22 


.0142 


.24 


2.76 
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TABLE V 



Results in stratified salt-water 

A 

Density gradient: .000020 #/ in 

Water depth: 19.5 inches 

Cylinder diameter: 1.5 inches 

Height of cylinder above bottom: 6.5 inches 



Point 


Force on 
cylinder 
(pounds) 


Amplitude 
of wave 
(inches) 


Inertia 

Coefficient 


1 


.0475 


.38 


2.64 


2 


.048 


.36 


2.78 


3 


.046 


.36 


2.68 


4 


.044 


.362 


2.53 


5 


.035 


.28 


2.60 


6 


.0322 


.25 


2.68 


7 


.029 


.24 


2.52 


8 


.0263 


.20 


2.74 


9 


.022 


.18 


2.54 


10 


.0198 


.16 


2.58 


11 


.075 


.61 


2.56 


12 


.072 


.58 


2.58 


13 


.068 


.54 


2.62 


14 


.064 


.50 


2.66 


15 


.0585 


.46 


2.65 


16 


.055 


.44 


2.60 


17 


.052 


.415 


2.61 


18 


.049 


.40 


2.56 


19 


.046 


.365 


2.62 


20 


.0365 


.29 


2.62 


21 


.028 


.22 


2.65 


22 


.023 


.18 


2.66 
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TABLE VI 



Results in stratified salt-water 

A 

Density gradient: .000020 #/in 

Water depth: 19.5 inches 

Cylinder diameter: 2 inches 

Height of cylinder above bottom: 6.5 inches 



Point 


Force on 
cylinder 
(pounds) 


Amplitude 
of wave 
(inches) 


Inertia 

Coefficient 


1 


.096 


.44 


2.56 


2 


,088 


.40 


2.35 


3 


.080 


.36 


2.60 


4 


.073 


.345 


2.48 


5 


.0645 


.29 


2.60 


6 


.0565 


.265 


2.50 


7 


.054 


.26 


2.44 


8 


.046 


.22 


2.45 


9 


.0392 


.18 


2.55 


10 


.036 


.16 


2.64 


11 


.0335 


.16 


2.46 


12 


.0295 


.135 


2.56 


13 


.0275 


.12 


2.68 


14 


.025 


.11 


2.66 


15 


.0225 


.10 


2.64 


16 


.146 


.66 


2.59 


17 


.107 


.47 


2.67 


18 


.078 


.35 


2.61 


19 


.0405 


.18 


2.63 


20 


.0265 


.12 


2.58 


21 


.094 


.42 


2.62 


22 


.038 


.17 


2.62 
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TABLE VII 



Results in stratified salt-water 

A 

Density gradient: .000042 #/ in 
Cylinder diameter: 1 inch 
Water depth: 19.5 inches 

Height of cylinder above bottom: 6.5 inches 



Point 


Force on 
cylinder 
(pounds) 


Amplitude 
of wave 
(inches).. 


Inertia 

Coefficient 


1 


.0134 


.21 


2.94 


2 


.0133 


.215 


2.85 


3 


.0104 


.18 


2.66 


4 


.010 


. 16 


2.88 


5 


.0093 


.15 


2.85 


6 


.008 


.13 


2.83 


7 


.0435 


.71 


2.82 


8 


.0405 


.67 


2.78 


9 


.0355 


.58 


2.82 


10 


.0265 


.44 


2.77 


11 


.026 


.42 


2. 85 


12 


.0245 


.40 


2.82 


13 


.0232 


.38 


2.81 


14 


.022 


.36 


2.81 


15 


.020 


.32 


2,88 


16 


.019 


.30 


2,91 


17 


.0172 


,28 


2,83 


18 


. 0154 


.25 


2,83 


19 


.014 


.23 


2,60 


20 


.0132 


.22 


2,76 
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TABLE VIII 



Results in stratified salt-water 

4 

Density gradient: .000042 #/in 
Cylinder diameter: 1.5 inches 
Water depth: 19.5 inches 

Height of cylinder above bottom: 6.5 inches 



Point Force on Amplitude Inertia 

cylinder of wave Coefficient 

(pounds) (inches) 



1 


.06 


2 


.058 


3 


.055 


4 


.052 


5 


.050 


6 


.046 


7 


.0435 


8 


.0417 


9 


.0403 


10 


.0375 


11 


.0356 


12 


.0325 


13 


.0295 


14 


.026 


15 


.025 


16 


.0158 


17 


.0108 


18 


.084 


19 


.063 


20 


.049 



.47 


2.60 


.44 


2.69 


.43 


2.61 


.40 


2.65 


.38 


2.69 


.36 


2.61 


.34 


2.61 


.33 


2.58 


.305 


2.69 


.295 


2.59 


.28 


2.59 


.26 


2.56 


.24 


2.51 


.21 


2.52 


.20 


2.55 


.13 


2.48 


.09 


2.45 


.63 


2.72 


.50 


2.57 


.39 


2.56 
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TABLE IX 



Results iii stratified salt~water 

Density gradient} * 000042 #/in 4 
Cylinder diameter} 2 inches 
Water depth: 19.5 inches 

Height of cylinder above bottom: 6.5 inches 



Point 


Force on 
cylinder 
(pounds) 


Amplitude 
Of wave 
(inches) 


Inertia 
Coeff ifcidftt 


1 


.141 


.64 


2*54 


2 


.130 


*58 


2.58 


3 


*110 


* 5 o 


2.53 


4 


.103 


*45 


2*63 


5 


.092 


*41 


2*58 


6 


.078 


*36 


2*49 


7 


.069 


*315 


2*52 


8 


.060 


.27 


2.56 


9 


.052 


*245 


2*44 


10 


.042 


*20 


2*42 


11 


.0385 


.18 


2*46 


12 


.036 


ill 


2*44 


13 


*0315 


*15 


2.42 


14 


.100 


*45 


2.56 


15 


.080 


*365 


2.52 


16 


.067 


.29 


2*66 


17 


.0565 


.26 


2*50 


18 


.0485 


*22 


2*535 


19 


*037 


.16 


2*66 


20 


*026 


*12 


2*49 
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TABLE X 



A comparison of theoretical and experimental inertia coefficients. 

Experimental values are means of numerous points taken during several 
runs. Note the very slight effect of the cylinder size and density 
stratification on the theoretical values of inertia coefficients. 







Inertia 


coefficients 


Cylinder 


size 


Theoretical 


Experimental 


(inches 


) 


value 

Fresh water 


mean 


1.0 




2.000 


1.998 


1.5 




2.000 


1.910 


2.0 




2.000 


2.094 




First 


4 

stratification, gradient 0.00002 #/in 


1.0 




2.500275 


2.78 


1.5 




2.500413 


2.61 


2.0 




2.500550 


2.565 




Second 


stratification, gradient 


0.000042 #/in 4 


1.0 




2.500715 


2.834 


1.5 




2.501074 


2.508 


2.0 




2.501433 


2.509 
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VI. Discussion of Results and Conclusions. 



The experimental data tabulated in tables I-X have shown that the 
inertia coefficient for circular cylinders depends on the density of the 
fluid and to a lesser extent on the density gradient and the diameter of 
the cylinder. For reasons yet unknown the inertia coefficients for the 
one inch cylinder are about 15% higher than those for the 1.5 and 2 inch 
cylinders. The values calculated with the appropriate values for the 
density, density gradient, and the diameter of the cylinder through the 
use of equation (24) gave an inertia coefficient of 2.5 as shown in 
Table X. This value varied very slightly with either the density gradi- 
ent or the cylinder size. It is, therefore, apparent from a comparison 
of the theoretical value with those obtained experimentally that the 
theoretical value is somewhat lower than those predicted experimentally. 
Among several reasons which could be provided, one is that the stratified 
fluid medium is confined and that this confinement, unlike in the case 
of homogeneous fluids, affects more adversely the blockage effect cited 
earlier. Additional errors may have come from the deviations of the 
density gradient from uniformity and from the presence of the connecting 
arms extending through a stratified media. Nevertheless, the agree- 
ment between the theory and experiments is sufficiently good to conclude 
that the inertia coefficient for a circular cylinder immersed in a strati- 
fied fluid media is about 25% larger than that in a homogeneous media. 

The consequences of the increase of the inertia coefficient are rather 
apparent when one considers the fact that in such media as the oceans, 
and in particular in layers close to the thermocline, the acceleration 
of a body requires a proportionally larger amount of power. 
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In the foregoing we have discussed so far the theoretical and 
experimental values of the inertia coefficients as they were calcu- 
lated or measured, but did not elaborate on the reasons as to why such 
an increase should occur. A careful analysis of the flow about a circu- 
lar cylinder immersed in a stratified fluid medium (See Fig. 6) shows 
that there are six instead of two stagnation points about the cylinder. 




Figure 6. Flow of a stratified fluid about a circular cylinder. 
Furthermore, there are on both sides of the cylinder two roughly tri- 
angular eddies symmetrically located with respect to either axis. It 
is apparent that the presence of additional stagnation points and the 
wedge shaped fluid body in the front and rear of the cylinder would 
necessarily increase the added mass and hence the inertia coefficient. 

In conclusion, it suffices to state that the inertia coefficient 
for cylinders immersed in a stratified fluid medium depends on the 
size of the cylinder, the density gradient, and the density of the 
fluid, and that the inertia coefficient is about 25% larger than that* 
of a cylinder immersed in an homogeneous medium. 



51 



VII. Recommendations for further Studies. 



It is strongly recommended that the project described herein be 
extended to flat plates, spheres, and other three dimensional bodies 
through the use of the same experimental apparatus. Even though an analy« 
sis of the stratified flow about a three dimensional body may be rather 
difficult, it is believed that an analysis similar to that presented 
here may be carried out for a flat plate vertically placed in the reser= 
voir, i.e,, normal to the direction of the horizontal component of the 
velocity. 
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